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ABSTRACT 

The economic feasibility of utilizing an ammonia bottoming cycle 
to improve the efficiency of a power plant is studied. The operating 
parameters of the bottoming cycle are examined to determine the optimum 
design binary cycle from a lifecycle costs viewpoint. 

The performance of the combined cycle is compared with that of a 
representative 'modern' steam plant: TVA's Bull Run. The analysis is 
based on current economic factors and the climatic conditions at the 
Bull Run site in Clinton, Tennessee. The optimum combined cycle is 
$98 million more profitable than the 'Bull Run' plant over the expected 
lifetime of the plant. The combined cycle also consumes 2 . 7 % less fuel 
than the unmodified steam plant. 

Optimization results are also presented for a range of environmental 
conditions at the plant location. 

The use of an ammonia bottoming cycle is economically as well as 
technically feasible provided the ambient water temperature is below 
69°F. 
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I. INTRODUCTION 



During the past few years, as the demand for electrical power has 
increased, mankind has become increasingly concerned over several 
important factors which will have impact on how electrical power will 
be generated in the future. The first of these considerations is the 
dramatic increase in the price of fuel and the realization that our 
fuel resources are limited. Secondly, we have become more concerned 
over how the electrical power generating plant affects our environment. 
The specific environmental considerations are how the fuel combustion 
process pollutes the air and how the heat which the plant rejects affects 
the temperature of the heat sink, usually a lake or river. 

In order to partially overcome these problems a plant with a higher 
thermal efficiency than presently exists is desired. The higher effic- 
iency plant would, for the same power output, burn less fuel and there- 
fore exhaust fewer combustion products to the atmosphere. Furthermore, 
as the plant efficiency increases the amount of heat rejected per unit 
power output decreases. 

The problem is: How to increase the efficiency of generating 

electric power? 

There are many ways to design a highly efficient power plant. Some 
of the proposed 'advanced' cycles are both expensive and exotic. The 
possibilities include: 

(1) Fuel -cel Is 

(2) MHD 

(3) Combined steam/gas turbine cycles 
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(4) Liquid metal topping cycles. 

(5) Bottoming cycles utilizing a refrigerant-type working fluid. 

In this thesis, the bottoming cycle will be examined. 

The thermal efficiency of an ideal cycle n^h " (Equation 1) 

where Tl is the average temperature at which heat is rejected and 
is the average temperature of heat addition to the cycle. 

The purpose of a bottoming cycle is to increase the thermal efficiency 
of a cycle by lowering T|_. Conventional steam plants are unable to take 
full advantage of a low heat rejection or condensing temperature due to: 

(1) The large specific volume of steam at low condensing 
temperatures and pressures. 

(2) The increased moisture content in the form of droplets 
as the condensing temperature drops. 

The above considerations require large, costly, and relatively 
inefficient low pressure turbine stages. 

In order to overcome these limitations and to take full advantage 
of a low temperature heat sink, another Rankine cycle using a different 
working fluid could be added to a steam plant. This 'bottoming' cycle 
would receive the energy which the steam plant rejects and would have 
the capability of lowering condensing temperatures. 

The selection of a working fluid for the 'bottoming' or sub-position 
cycle cannot be accomplished by a closed form equation, but rather by a 
series of trade-offs. The desired properties of the working fluid are; 

(1) Have a relatively low specific volume at the condensing 
temperatures to be considered. 

(2) Be a liquid at ambient temperatures, so it can be pumped. 

(3) Be non- toxic and non-corrosive. 
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(4) Have a low specific heat and a large latent heat of 
vaporization over the temperature range to be con- 
sidered. (This results in the average temperature of 
heat addition being as high as possible). 

(5) Have a condensing pressure above atmospheric pressure. 

(6) Have a high speed of sound at turbine exit conditions. 

(7) Have good heat transfer properties. 

(8) Be chemically stable and non-flammable. 

(9) Be inexpensive and plentiful.^ 

No known substance satisfied all of these requirements, but 
Ammonia (NH3) seems to be the best overall choice. 

The feasibility of the ammonia bottoming cycle will be evaluated 
by investigating a representative 'modern' steam plant modified to 
include the sub-position cycle. The thermal efficiency, capital, 
and operating costs of the binary cycle plant will be compared with 
those of the unmodified 'modern' plant. With these factors in mind 
an optimum design combined cycle plant will be generated, which will 
result in maximum savings to the owner/operator over the expected 
lifetime of the plant. 
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II. LITERATURE SURVEY 

2 

The concept of a bottoming cycle is not new. In 1961, Aronson 
proposed such a binary-vapor cycle, utilizing Freon -12 as the working 
fluid. His analysis showed a heat rate improvement of 1-2% over that 
of a 'modern* steam plant as well as establishing the economic feasi- 
bility of the cycle. 

3 

In 1969, Wood stated that a binary-vapor bottoming cycle was not 
feasible, citing available energy loss due the temperature difference 
across the steam condenser/secondary fluid boiler. 

The Aronson proposal has been supported by Steigelman et al^, who 
demonstrated the desirability of a binary-cycle using cooling towers 
provided the ambient air temperature was low enough. Furthermore, El 

5 

Ramly and Budenholzer professed the advantages of a bottoming cycle to 
operate in conjunction with a nuclear power plant. 

Perhaps the most comprehensive study to date has been that of 
Slusarek®, in which the feasibility of the ammonia bottoming cycle 
was studied in detail. The Slusarek report included a thorough component 
design, an economic analysis, and an engineering optimization. While 
the majority of the report is accurate, two significant shortcomings 
must be noted: 

(1) The heat transfer coefficients in the ammonia boiler were 
calculated using inappropriate correlations. 

(2) The impact of the steam condensing temperature and the 
temperature difference between the main and bottoming cycles 
were not adequately examined. 
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In this study, the operating parameters of a bottoming cycle will 
be studied and optimized in order to prepare a preliminary design of the 
best plant from an engineering as well as economic viewpoint. 



III. THE COMBINED CYCLE 

As outlined in chapter one, the bottoming cycle utilizes the heat 
rejected by a steam plant as its energy source. A simplified diagram of 
a combined steam-ammonia plant, hereafter called the 'combined' cycle, 
is shown in figure one. 

To carry out the analysis of the feasibility and economic optim- 
ization of the combined cycle, an existing 'modern' steam plant was modified 
to include a bottoming cycle. The Tennessee Valley Authority's 'Bull 
Run' plant was selected for this role. 

Bull Run is a 900 MW, coal fired, supercritical steam plant which 
first went into operation in the mid 1960's. The plant features a 
boiler exit temperature of 1000°F and a single reheat to the same tem- 
perature. The design condenser pressure is 1.5" of mercury, which 
corresponds to 91.7°F and .74 psi. The thermal efficiency of Bull 
Run, excluding the boiler losses is 45.78%. Including the boiler 
efficiency of 89%, the cycle efficiency is 40.74%.^ 

The Bull Run plant flow diagram and the Temperature vs. Entropy 
diagram are shown as figures two and three respectively. 

The ammonia bottoming cycle is a simple Rankine cycle. Feedheating 
must be utilized since the feasibility analysis is sensitive to small changes 

O 

in the efficiency of the bottoming cycle. 

Based on earlier studies, a four stage turbine for the ammonia cycle 
seems to be a logical choice. Consequently, the cycle will have three stages 
of feedheating. Data from steam plants indicates that the maximum efficiency 



- 17 - 



o 

CVJ 




>- 

LU 




- 18 - 



FIGURE 1: SIMPLIFIED DIAGRAM OF THE COMBINED CYCLE AMMONIA 
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ENTROPY 



FIGURE 3 TEMPERATURE - ENTROPY DIAGRAM OF TVA'S BULL RUN PLANT 
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STEAM NUMBERS REFER TO FIGURES 2 and 3 
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INCLUDES FEEDHEATING 



increase for three stage feedheating occurs when the working fluid 
temperature is elevated three-quarters of the difference between the 

9 

condensing and evaporating temperature by feedheating. 

The Temperature-Entropy diagram of the bottoming cycle is shown 
as figure 4. Note that the ammonia is not superheated. In order to 
maximize the efficiency of the combined cycle, the average temperature 
of heat addition to the ammonia cycle must be as close as practical to 
the steam condensing temperature. If superheat were introduced, the 
average temperature of heat addition for the bottoming cycle would 
drop, negating the efficiency gains usually derived from superheating. 

In order to bring the steam and ammonia cycles together into one 
combined cycle certain modifications must be made to the two cycles. 

The modifications are: 

(1) Replacing the water-cooled steam condenser with a steam 
condenser/ammonia boiler. 

(2) , Raising the backpressure of the low pressure steam turbine 

and combining the remaining LP turbine stages with the IP 
turbine. 

(3) Removal of feedheating stages in the steam plant made un- 
necessary by raising the steam condensing temperature. 

The T-S diagram and component flow chart for the combined cycle 
are shown as figures five and six respectively. 

Since the high temperature end of the steam plant is not altered, 
the only independent parameter in the steam side of the combined plant 
is the steam condensing temperature, T^. On the ammonia side there are two 
independent variables: 

(1) The ammonia boiling temperature TF(1). 

(2) The ammonia condensing temperature, TF(5). 
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ENTROPY 

FIGURE 4 T-S DIAGRAM OF BOTTOMING CYCLE 



STATE 

a 

b 



c 



d 

FLOW 
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DESCRIPTION 

Condenser exit (saturated liquid) 

Exit of third feedheater, boiler inlet (liquid) 
Boiler exit, turbine inlet (saturated vapor) 
Turbine exit, condenser inlet (liquid and vapor) 
DESCRIPTION 

Mass flow rate through boiler 

Mass flow rate through condenser 

Mass flow rate of feedheating extraction #1 

Mass flow rate of feedheating extraction #2 

Mass flow rate of feedheating extraction #3 



TABLE 2 DESCRIPTION OF T-S DIAGRAM (FIGURE 4) 
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FIGURE 5 T-S DIAGRAM OF COMBINED CYCLE 
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FIGURE 6 COMPONENT FLOW CHART FOR COMBINED CYCLE 




TABLE 3: COMBINED CYCLE DESCRIPTION 



STREAM NUMBERS REFER TO 

STREAM 

Fuel 

Air 

Exhaust 

Water 

Supercritical Liquid 

Steam 

Steam 

Steam 

Steam/Water 

Ammonia 

Ammonia 

Ammonia 

Ammonia 

Water 

Water 

* INCLUDES FEEDHEATING 



FIGURES 5 AND 6 

NO. 

a 

b 

c 

1 

2 * 

3 

4 

5 

6 
7 

8 * 

9 

10 

11 

12 



DESCRIPTION 
Fuel inlet 

Air inlet for combustion 
Combustion products 
Condenser exit 
Boiler inlet 

Boiler exit/HP turbine inlet 
HP turbine exit/reheater inlet 
Reheater exit/ IP turbine inlet 
IP turbine exit/condenser inlet 
Condenser exit 
Boiler inlet 

Boiler exit/NH 3 turbine inlet 
NHj turbine exit/condenser inlet 
Condenser coolant water inlet 
Condenser coolant water exit 
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The three variables mentioned above are the only independent 
variables in the analysis. The significant dependent variables are: 

(1) The entropy averaged temperature of heat addition to the 
ammonia cycle, TAV, which is a function of TF(1) and TF(5). 

(2) The difference between the steam condensing temperature T-, 
and TAV. This variable is called Aj and is a function of^ 
all three independent variables. The magnitude of this 
parameter is related to the available energy lost due to 
the temperature drop between the cycles. 

(3) The Log Mean Temperature Difference (LMTD) in the steam 
condenser/ammonia boiler. This factor which is numerically 
close to Af, is a function of the same three variables, 
and is a significant consideration in the sizing and cost 
estimation of the steam condenser/ammonia boiler. 

The object of the combined cycle optimization is to consider 
how the three independent variables affect the plant performance and 
cost. 
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IV. OPTIMIZATION PRELIMINARIES 

In this chapter the impact of the independent and dependent 
variables on the analysis will be examined. In addition, limita- 
tions on the range of these parameters will be discussed where 
appropriate. 

a. The ammonia condenser and the ammonia condensing temperature. 

In order to maximize the overall cycle efficiency, the lowest 
possible ammonia condensing temperature, TF(5), is desired. However, 
for a power plant with a water-cooled condenser there are three factors 
which limit how low the condensing temperature can be. The factors 
are: 

(1) The condenser coolant water temperature at inlet must be 
greater than 32°F. 

(2) The water velocity through the condenser must not be so 
high as to cause excessive corrosion in the tubes. 

(3) The amount of power required to pump the coolant water 
through the condenser tubes must not negate the advantages 
of lowering the condensing temperature. 

The reason why the water temperature at inlet must be greater 
than 32°F is that it must be a liquid in order to be pumped. If the 
ambient temperature is below 32°F, a cooling tower would be utilized 
in place of a wet condenser. The considerations of water velocity and 
pumping power are not quite so obvious, but are related to each other. 

In a heat exchanger, such as a condenser, the temperature difference 
which will partially determine the condenser size is 
AT = Tc - {Tj + Tq)/ 2 (equation 2) 

where Tq is the condensing temperature, Tj is the; water inlet temperature. 
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and Tq is the water exit temperature. The temperature vs. length plot 

for a condenser is shown figure 7. 

The heat transferred in the condenser Q may be expressed as: 

Q = UA AT (equation 3) 

where A is the heat transfer area and U is the overall heat transfer 

coefficient for the heat exchanger. The overall heat transfer coefficient 

may be computed by combining the thermal resistances in the unit as if 

they were resistors in a parallel electrical circuit. Assuming the 

thermal resistance of the condenser tubes to be negligible, and the 

thermal resistances for convection and condensation to be the 

inverse of the average heat transfer coefficients on their respective 

sides of the unit, then U = l/(l/bQondensing ^ ^/^convection on) (equation 4) 

fluid water side 

For shorthand purposes the following symbolism will be used: 

^condensing fluid ->■ *^NH 3 cond f^convection on water side f^water 
In terms of the fluid and thermal parameters for a horizontal tube 
condenser: 

hNHi = .728 [1+.2C AT(n-l)/hfg] ?tg(p)(p " Pv)k3 h'fa)/(nDy AT„) (equation 5) 
cond ^ 

^water ' -^23 kf/D (pfVD/yf)*® (Cpp/kf)-^ (equation 6) 

where C = specific heat of condensing liquid 

Cp = specific heat of water 

D = tube diameter 

g = acceleration of gravity 

hfg = heat of vaporization of condensing fluid at condensing 

pressure 

h'fg = l+(.68 C AT^)/ hfg 

k = thermal conductivity of condensing liquid 

k^ = thermal conductivity of water 
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LENGTH FROM WATER INLET 

FIGURE 7 TEMPERATURE VS. LENGTH FOR A CONDENSER 
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n = number of tubes in vertical row 

p = density of condensing liquid 

Pf = density of water 

Py = density of condensing vapor 

y = dynanic viscosity of condensing liquid 

Pf = dynamic viscosity of water 

V = water velocity through the tubes 

AT = temperature difference between condensing fluid and 
tube wall 

Furthermore, an energy balance for the condenser states: 

Q ^ water Swater ^^0 ‘ = *"condensing fluid ^e ^fg (equation 7) 

and solving for Tq - Tj, the water temperature increase: 

^0 “ condensing fluid ^e *^fg)^(”V/ater ^Pjyater^ (equation 8) 

where Xg is the condensing fluid quality. Combining equations (2), (7), 
and (8) then, 

at = Tq - (Tj + Tq) /2 = Tc - (Tj + Q/2my,afgy.) (equation 9) 

From equation (3) it may be seen that for a condenser to be of 
minimum size and cost once the heat transfer amount Q has been selected, 
the product UAT must be as great as possible. 

However, from equations (5) and (6), it may be shown that AT changes 
much more rapidly than U, apparently stating that AT should be maximized. 
Equation (9) indicates that once Tq is chosen (and ideally as low as possible), 
there are only two ways to increase AT. The first method is to decrease 
Tj, the water inlet temperature, which is difficult since Tj is a function 
of the plant site. 

On the other hand, if the Q/2my,afgj. term in equation (9) is small 
the value of AT will remain large. The problem is that the only way to 
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decrease the Q/2ir\^ater ^ to increase the mass 

flow rate of water through the condenser tubes. This brings one 
back to the problems of water velocity and pumping power stated earlier 
in this chapter. 

In order to keep the economic analysis of the components as 
simple as possible, the combined cycle plant condenser under consider- 
ation will have the same ratio of cycle power output to surface area 
as the 'Bull Run' unit. This assumption is valid, since the condenser 
limitations are on the water side. The 'Bull Run' condensing system 
consists of 41368 tubes, each of which has an outer diameter of .875", 
and a total heat transfer surface area of 320,000 feet . Since the 



(where Af|ow perpendicular to the direction of water flow, 

n is the number of tubes, and D is the tube inner diameter), then 



Now that an expression for the water velocity as a function of 
the mass flow rate and the condenser geometry has been defined, a limit 
on the mass flow rate and TF(5) can be specified once a limitation on 
the water velocity is established. 

Since ammonia reacts readily with copper, stainless steel is the 
likely candidate for the tube material. Studies indicate that for steel 
tubes corrosion will not be excessive as long as the water velocity is 
below 15 ft/sec.^^ 



mass flow rate, Pwater ^mter ^flow 



(equation 10) 
(equation 11 ) 




'^water " %ater/Pwater ^flow ~ Aifiwater/Pwater nTrO (equation 12) 
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This water velocity may require excessive pumping power which 
would be counter-productive. The pumping power, PP may be expressed 
as: PP = (AP X Q)/rip (equation 13) 

where Q is the volumetric flow rate, Op the pump efficiency, and 
AP is the pressure drop through the tubes. 

The volumetric flow rate Q, is 

^ water/p water (equation 14) 

and the pressure drop, AP is: 

AP - 4 f (L/D) (py^g'^gP V vyater/^^o) (equation 15) 

where f is the pipe friction factor, L the length and g^ is the 
units correction factor. 

For turbulent flow in tubes: 

f = .079/Rg'^^ (equation 16) 

where is the Reynolds number and is equal to 

"e = ((-water t»/ "water (equation 17) 

Combining equations (13) through (17), the total pumping power 

is: 

PP - .316 / (o„V„D/U„)-26 (L/0)(0^V^ / 2 g^) (l"„/p„)(l/np) (eq. 18) 
recalling that L = TOTAL HEAT TRANSFER AREA (A)/(n7rD) (eq. 19) 
and m = mrD‘^)/4 (equation 20) 

the pumping power expression simplifies to: 

PP = (.079 A V^*^^)/ (4 9^ Op) (equation 21) 

in which V is the most significant factor and the only one that will 
change noticeably from one set of parameters to another. 
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Thus far it has been shown how the ambient water temperature 
Tj, and practical engineering considerations limit the ammonia 
condensing temperature. 

To summarize: 

(1) Once a plant site has been chosen 
Tj is fixed. 

(2) For a fixed Tt and a given condenser, 
the mass flow^rate (and velocity) must 

be maximized in order to keep the ammonia 
condensing temperature as low as possible, 

(3) Limitations on the water velocity are 
imposed by corrosion considerations. 

(4) In addition to #3, the water velocity must 

not be so large as to require excessive pumping 
power which would negate the advantages of 
lowering TF(5), the NH^ condensing temperature. 

The interface between the gains of lowering TF(5) and the 
pumping power requirements will be considered later. 

b. Steam Condensing Temperature, T^ 

It is desired to find the steam condensing temperature which 
will give the maximim cycle efficiency. To obtain the optimum T^, 
it must be determined whether the steam or ammonia cycle makes better 
use of the low temperature end of the cycle. 

Assuming that the steam and anmonia turbines have the same 
isentropic efficiencies in the two-phase region, the deciding factor 
in the selection of T^, is the power loss due to moisture droplets in 
the turbine. 

Common practice indicates that the power loss due to wetness is 
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one percent for each percent of average wetness in that turbine stage. 

Applying this wetness consideration to the 'Bull Run' plant, a 
plot of power loss due to wetness versus the steam condensing temperature 
may be created as shown in figure 8 . 

Once the ammonia condensing temperature TF(5) has been selected, 
the power loss analysis may be carried out for the ammonia cycle as 
well. However, since the third independent variable, the NH^ boiling 
temperature TF(1), has not been fixed the wetness calculation must be 
carried out for each value of TF(1). This procedure produces a family 
of curves one of which is shown in figure 9. 

Figure 8 demonstrates that the power loss for the steam cycle 
increases as T 5 drops, and figure 9 shows the opposite trend for the 
ammonia cycle. 

The optimum value of T^ may be found by superimposing figures 8 
and 9 and plotting the combined cycle wetness loss versus T 5 . 

Figure 10 shows the superimposition of figures 8 and 9 and the 
resulting total wetness loss. Note that a mimimum loss occurs around 

O 

T 5 = 131 F for the parameters chosen. Furthermore, if the combined 
cycle efficiency is plottedversus the non-dimensionalized temperature 
difference between the steam and ammonia cycle for various values of 

O 

Tg, as shown in figure 11, it is clear that T^= 131 F is the optimum 
choice. 

The selection of a steam condensing temperature at this level is 
different from the studies cited in chapter II. All of the earlier work 
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assumed a value of in the neighborhood of 220° - 250°F without 
considering the wetness losses. 

The analysis which generated figures 10 and 11 will be discussed 
in the next chapter. 

c. The Ammonia Boiling Temperature TF(1) and the Steam Condenser/ 
Ammonia Boiler. 

The third independent variable is TF(1), the ammonia boiling tem- 
perature. This variable is most significant in its impact on the cycle 
efficiency and the size and cost of the steam condenser/ ammonia boiler. 
The value of TF(1) is not as important as its relationship with the steam 
condensing temperature T^. This relationship defines the dependent 
variable AT, the average temperature difference between the steam and 
ammonia cycles. 

As indicated in chapter three and graphically shown in figure 11, 
the combined cycle efficiency is highest when AT is as small as possible. 
Unfortunately, since the average temperature difference AT is numerically 
similar to the log mean temperature difference in the steam condenser/ 
ammonia boiler and the heat transferred in the heat exchanger Q is: 

mean- (equation 22) 

Thus as the temperature difference decreases the heat exchanger area 
increases. 

It is obvious that the efficiency gains derived by a small AT may 
be negated by the additional capital cost of a large steam condenser/ 
ammonia boiler. 
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The task at hand is twofold: 

(1) Determine how aT affects the overall thermal efficiency. 

(2) Determine how AT affects the steam condenser/ammonia 
boiler size and cost. 

Task one will be considered in the next chapter. 

In order to determine how AT affects the heat exchanger size and 

cost, a detailed heat transfer analysis of the unit must be made. 

% 

The steam condenser/ammonia boiler to be considered is a shell 
and tube arrangement with the steam condensing outside the array of 
horizontal tubes through which the ammonia flows. 

For analytical purposes the SC/AB may be divided into two zones. 

In the first zone, the ammonia is heated from the feedheating outlet 
temperature, TF(2) to the ammonia boiling temperature TF(1). The heat 
transfer mechanism on the ammonia side in this region is forced con- 
vection, with subcooled boiling effects being neglected, and on the 
steam side the mechanism is condensation. 

In the second zone, the ammonia temperature remains constant as 
the ammonia is heated from a saturated liquid to a saturated vapor state. 
The heat transfer mechanism on the NH^ side in this zone is very complex, 
and includes effects of nucleate boiling and forced convection. The heat 
transfer mechanism on the steam side is once again condensation. 

In both the non-boiling (zone 1) and boiling (zone 2) regions the 
heat transfer coefficient on the steam side hg^gg^ may be expressed as: 

^steam ' [1 + .2C AT(n-l)/hfg] .^g(p)(p- p^)k^ h'fg)/nDu AT^) (eq. 23) 

where the equation and symbols are identical to equation 5, which was used 
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in the NH^ condenser calculations. 

On the ammonia side in the non-boiling region the forced convection 

8 4 

coefficient is: hf^^^ = .023 (kf/D) (p^VD/y)* (Cpy/k^)' (eq. 24) 

which is the same expression as equation( 6), except that the fluid 

properties for ammonia are used here. 

In zone 2, the boiling region, the heat transfer coefficient may 

be calculated by using several empirical relationships. The most 

12 

reliable seems to be the Chen correlation which states that the 
heat transfer coefficient in the two-phase region where boiling and 
forced convection are present, may be calculated as follows: 



^^TWO PHASE = I^NUCLEATE BOILING ^ ^ + ^poRCED CONVECTION 2 PHASE 25) 

where S is the boiling suppression factor, a measure of how much 

the fluid motion upsets the normal boiling mechanism and F is the 
Chen factor which is a measure of how much the fluid motion induced by 
boiling upsets the normal forced convection mechanism. 

The Chen factor may be calculated as follows: 

(1) Determine the Martinelli -Nelson factor X^-^, where 

^tt ~ (Pf/Pg) ((l-x)/x) (equation 26) 



Note that p and p represent the density of saturated vapor 
and liquid respectively, y^ and y the viscosity of saturated 
liquid and vapor respectively, ana x is the thermodynamic 
quality. 



(2) Utilizing figure 13, which displays F versus read off F. 

The boiling suppression factor S may be found as follows: 

(1) Find the two-phase Reynolds number which is represented by 
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(3) Utilizing steps one and two the forced convection of ammonia 
heat transfer coefficient may be found. 

(4) Solve the heat flux equation by 

Qnb/% = h steam (Ts-Twall) = hNH 3 FC ('•'wal 1 -TNH 3 ) (eq. 33) 

iteration. The reason that this relationship cannot be solved 
in closed form is that the expression for hg^gg^^ (equation 23) 

contains the term (T$-T,^g^]) ^ making a numerical result 

preferable. The iteration is carried out by assuming a tube 
wall temperature and incrementally changing T^^g]] until h^^ggj^ 

(^S-Twall) equals •^NH 3 p^^wall'^NH 3 ) • *^°te that T|y|H 3 is the 

numerical average between the ammonia inlet and outlet temperature. 

(5) Divide the non-boiling heat transfer by the heat flux to 
determine the area of the non-boiling section, A|\^g. 

( 6 ) Calculate the pressure drop in this section of the SC/AB and 
determine that the design choices made in step one are 
responsible. Use equation 15 and substitute ammonia properties 
for those of water. 

The procedure for the boiling section calculations are basically the 
same, but are more complex. 

Note that in equation(26)and(30) ,the variable 'x' which represents 
the thermodynamic quality appears. In order to calculate the heat 
transfer coefficient for a region where the quality goes from zero to one, 
what value of x should be used? A simple average between the end values 
is not suitable since in both equations(26)and(30) ,the 'x' term is non- 
linear. However, in order to make this problem less difficult to solve, 
an approximation which yields good accuracy is to use an average quality 
in the analysis, provided the change of the quality is small. 

In keeping with this idea, the boiling region of the SC/AB is 
divided into ten segments, each having a quality change of one-tenth and 
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and utilizing the average quality in that segment in equations (26) and 
(30). 

The configuration of the steam condenser/ammonia boiler is determined 
(step one for non-boiling section) and the amount of heat transfer in 
the boiling region, Q^, as well as the NH^ mass flow rate (step 2) 
are fixed. Then for each quality segment of the boiling region the heat 
flux iteration must be performed: 



''b segment 




S + ^'2PHASE 

FORC. CONV 



n(TF(l)-Twaii) (eq. 34) 



This iteration is more complex than the non-boiling case since hNucLEATE 

BOILING 

and both contain temperature difference terms raised to fractional 

powers . 

Once the heat flux for each segment has been determined, the heat 
transfer area of that segment may be found. The total boiling region area 
A is the total of the areas of the ten quality segments. 

D 

Finally, the pressure drop in the boiling region must be evaluated 
to determine design feasibility. Using a homogeneous flow model^^, the 



pressure drop for the boiling region is: 

AP = 2fjp pv\ ^ o\l^ (Vfg/Vf) 

D 2vf 

where fjp = .079/ (pVD/p)*^^ 



(equation 35) 



(equation 36) 



and p is the average viscosity and may be expressed by the McAdams cor- 



relation^^: p- = xp + (l-x)p . 

9 f 



(equation 37) 
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Once again quality is important and the calculation of the pressure 
drop is carried out in the same manner as the heat transfer coefficients, 
by assuming the calculation to be piecewise linear, 
d. Chapter summary 

In this chapter, the impact of the three independent variables 

TF(5), TF(1), and on the heat exhanger design and cycle efficiency 
have been examined. Specifically, the NH^ condenser performance and 
pumping requirements, the selection of an optimum T^, and the sizing 
of the steam condenser/ammonia boiler have been considered. In the next 
chapter, the thermodynamic analysis will be discussed and the impact of the 
three variables on cycle efficiency. 
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V. THERMODYNAMIC ANALYSIS OF THE COMBINED CYCLE 

In this chapter, the impact of the three design parameters, the 
ammonia boiling temperature, ammonia condensing temperature, and steam 
condensing temperatures on the combined cycle efficiency will be studied, 
a. The Steam Cycle 

The steam condensing temperature, T^, is the only variable which 
affects the thermal efficiency of the steam cycle. Specifically, it 
would be desirable to know how the steam cycle efficiency ng-p, varies 
with Tg, the steam condensing temperature. 

The power output of the steam cycle is affected by two considerations 
as the steam condensing temperature is increased to accomodate the bottoming 
cycle. The first of these considerations is that the LP steam turbine 
produces less power due to the higher exhaust backpressure. Secondly, 
since the steam is being condensed at a higher temperature than the 
'Bull Run' plant, less steam must be extracted from the power producing 
turbines for feedheating purposes. Note that those considerations 
affect the cycle efficiency in opposing manners. 

With regard to the loss of power due to increased backpressure, 
only the low pressure steam turbine is affected. The power loss associated 
with the higher backpressure is: PL = liisteam ( l^i - hg ) (equation 38) 

where h^- is the enthalpy of the stage exhaust steam for the increased back- 
pressure case and hg is the stage exhaust enthalpy for the unmodified 
Bull Run plant. Equation (38) must be modified to include the difference 
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in wetness losses between the original and the modified steam plant. 

? 

Accordingly if.h^. is theinlet enthalpy, x, is the average quality for 
the modified plant, and Xq is the average quality for the base case, J 
equation (38) becomes: f 

PL - m steam ( (^^i " *^o) ^o “ ^1^ ^l) (equation 39) ^ 

Equation (39) is valid for each stage in the LP turbine provided , 
that the proper steam flow rate, m steam used. 

In addition to the difference in power output due to the higher 
backpressure, the changes in the feedheating extractions also affect the 
power output. 

In the 'Bull Run' plant, the extracted steam is condensed in a 
closed-type feedheater, transferring heat to the feedwater until the* 
condensed extraction temperature is 5°F higher than the feedwater outlet 
temperature. Ultimately, the condensed extraction is mixed with the 

f ' 

feedwater from the main steam condenser. 

In order to determine the cycle power variation due to changes in 
feedheating extractions, an expression for the mass flow rate of steam 
'a' in figure 15 must be established. 

Referring to figure 15, an energy balance of the closed feedheater 
indicates: 

•"a (^a “ ^b^ ~ ^d (^e ~ ^d^ (equation 40) » 

A design constraint for the 'Bull Run' feedheaters states: 

Tjj = Tg + 5’’F (equation 41) 
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and since = 1 BTU/lbm - °R, then: 



(equation 42) 



hb = hg + 5 BTU/lbm 
An energy and mass balance in the mixer states: 
itid = ma + me (equation 43) 

l^d ^ ("^a *^b mchc)/(ma + m^) (equation 44) 

Combining equations (40) through (44), the extraction flow, mg, 

may be expressed as: 



ma mg 



(ma + >^)he - ma(he + 5) -rtichc 



(itig + mr)(ha - 2hp - 5) + ma(he + 5) + 



m, 



c c 



(equation 45) 



The condenser flow rate is also affected by the amount of 
feedheating extraction. The definitive relationship is: 



"*c = "iBOILER - "’OTHER EXTRACTIONS 



- m^ 



renaming meoRER “ "’OTHER EXTRACTIONS = "’g 
and combining equation (45), (46), and (47): 



(equation 46) 
(equation 47) 



™a ' <™g - "a' 



n>ghe ■ "aC'e * 5) - (n.g 



"’a)'’c 






(eq. 48) 



The variables in equation (48), m and h are fixed at the values 
of the unmodified plant. 

The feedheater exit enthalpy will be fixed as a design constraint 
in this analysis. Thus, there are only two unspecified variables, the 
first is the enthalpy of the main condenser condensate, hg> which is a 
direct function of the steam condensing temperature, T^. The second 
variable is the extraction mass flow rate, ihg. Based on the arguments 
presented, it is now possible to state that: mg = f(T^) (eq. 49) 
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where the function f is specified in equation (48) and in the 
relationship between h^ and T^. which may be found in the steam 
tables. 

Therefore, if the extraction mass flow rate for the unmodified 

cycle is liia , and the extraction flow rate for the new steam con- 
0 

densing temperature is irig^ , the power loss due to the higher con- 
densing temperature (including the effects of equation (39)) is: 

PL = - Iio)Xq “ (equation 50) 

where m is the turbine stage mass flow rate for the unmodified 
^0 

case, and m is the modified cycle mass flow rate. The modified 
^1 

cycle flow rate may be expressed as: 

>"si = msg + ('"aQ - iTig^) (equation 51) 

Combining equations (50) and (51), the power loss is : 

PL = ms^(h. - ho) Xo - (ms^j - (% " ~ ) x, (eq. 52) 

The thermal efficiency of the modified steam cycle becomes: 

nth =(^o ■ (equation 53) 

steam ”steam 

where Pq is the unmodified cycle power output and QA^^eam 
added to the steam cycle. 



Consequently, the heat rejected by the steam cycle, which is 
also the heat added to the bottoming cycle is: 

\team " ^83 " ^f^steam “ "^th^^gg^) (equation 54) 

It may be noted that effect of only one stage of feedheating was 
shown in the modification of the 'Bull Run' plant. For the range of 
values where the steam condensing temperature T^, yields optimum combined 
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cycle efficiency, only the first feedheating stage is affected. 

However, a more complex analysis was undertaken for values of 
which affected the second and third feedheating stages, 
b. The Ammonia Bottoming cycle 

The ammonia bottoming cycle is a Rankine cycle with three stages 
of feedheating and no superheating. 

In carrying out the thermodynamic analysis of the bottoming cycle, 
the following symbols are used: 
s = entropy (BTU/lbm °R) 
h = enthalpy (BTU/lbm) 

P = Pressure (psi) 

TF = temperature (®F) 

V = specific volume (ft /Ibm) 

f = used as a subscript, indicates saturated liquid state 

g = used as a subscript, indicates saturated vapor state 

fg = used as a subscript, indicates property change associated with 

vaporization or condensation 

Figure 16 is the T-S diagram for the bottoming cycle with the no- 
menclature to be used in the analysis. Note that since there are 
three stages of feedheating, there are five temperature levels in the cycle. 
The temperature levels are numbered sequentially, starting with level one 
which is the ammonia boiling temperature and culminating with level five 
which is the ammonia condensing temperature. 

The power extracted from each turbine stage is the product of the 
mass flow rate and the enthalpy drop in the stage, minus mechanical, kinetic 
energy and moisture losses. 
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FIGURE 16 TEMPERATURE-ENTROPY DIAGRAM WITH ANALYSIS 
NOMENCLATURE FOR NH 3 CYCLE 
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The ideal enthalpy drop is found by considering the expansion 
in the turbine stage to be isentropic. Therefore, using the first 
stage as an example: 

s(2s) = s(l) (equation 55) 

and the quality at state 2s is: 

x(2s) =(s(2s) - Sf(2))/s^g(2) (equation 56) 

therefore, the enthalpy of state 2s is: 

h(2s) = hf(2) + hfg(2) • x(2s) (equation 57) 

However, since in reality the expansion is not isentropic, nis> the 
isentropic efficiency is introduced: 

=(h(l) - h(2))/(h(l) - h(2s)) (equation 58) 

therefore the actual enthalpy drop in the first stage is: 

Ah]2 = h(l) - h(2) = n-jg (h(l) - h(2s)) (equation 59) 

The loss due to wetness corresponds to one percent power loss 
for each one percent of average moisture in that stage. Conversely, 
the wetness losses could be accounted for by multiplying the enthalpy 

drop by the average quality in the stage. For the first stage, the exit 

quality x(2) is: 

x(2) =(h(2) - h^(2))/ hfg(2) (equation 60) 

and the average quality in the first stage 

x-| = (x(l) + x(2))/2 (equation 61) 

The mechanical, stage, and exit losses are each assumed to be one 
percent. Expressing these three factors in one expression called 
mechanical efficiency, which is simply one minus the three one 
percertt loss, thus = .97. (equation 62) 
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Therefore, the power extracted from the first stage of the ammonia 
turbine Pj^ , is: 

PTi = mi(h(l) - h(2)) • • (x(l) + x(2))/2) (equation 63) 

The same analysis may be applied to the final three stages of the 
ammonia turbine, provided the proper mass flow rate for each stage is 
used. Recall that each successive stage has a smaller flow rate due 
to feedheating extractions. 

The purpose of feedheating is to improve the cycle efficiency by 
increasing the average temperature of heat addition to the cycle. 

This is accomplished by extracting some of the working fluid and 
utilizing it to heat the boiler feed. The type of feedheater used in the 
analysis is an open-type, in which the feed is mixed with the extracted 
fluid. 

The extracted fluid mixes with the feed after the feed pump exit. 

The mass flow rate of the extraction is adjusted so that the mixture is 
at the saturated liquid stage. For example, the condensate (state F5 of 
figure 16) is pumped to the pressure corresponding to saturation at tem- 
perature level four. Ammonia is extracted at state 4 and mixed with the 
feed until the mixture is at the state F4. This process is repreated until 
the feed is at the boiler inlet condition. 

An energy balance in each feedheater states: 

"^feed ^^feed “ *^extraction ^*^extraction (equation 64) 

’^extraction infeed Ahfgg^j/ (equation 65) 

By substituting, using the symbols from figure 16 equation (65) becomes 
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nig^ = (1 - me^) (hf(2) - (hf(3) + Wp(3))/(h(2) - hf(2)) (equation 66) 
iiie2 = (1 - - me2)(M3) - hf(4) + Wp(4))/(h(3) - hf(3)) (eq. 67) 

liies = 0 - rnei ' >^e2 ' ^ Wp(5))/(h(4) - hf(4)) (eq. 68) 

where irii = 1 and % = 1 - m. - ~ • 

1 ^ e-| e2 

The term Wp which appears in these three equations is the enthalpy 
rise associated with the pump work. The subscript number indicates 
the temperature level at the pump inlet. The pump work, assuming the 
specific volume change across the pump is negligible is: 

Wp (n) = v^(n) (P^^.| - P^)/hp X 144/778 BTU/lbm (equation 69) 

where n is the temperature level at the pump inlet, and np is the pump 
isentropic efficiency. 

The heat added to the ammonia cycle is the product of the mass 
flow rate and the enthalpy rise in the ammonia boiler. The heat added 
Qa NH 3 = - (hf(2) + Wp(2)) (equation 70) 

to the ammonia cycle is also the heat rejected by the steam cycle. 
Combining equations (54) and (70), 

NH3 = Qa steamd-nth steam^ = + Wp(2)). (eq. 71) 

The left hand side of this equation is fixed by the steam condensing 
temperature T^. The enthalpies on the right side are functions of the 
ammonia boiling and condensing temperatures, TF(1) and TF(5) respectively. 
Therefore, once these three variables are fixed, iri] the mass flow rate 
through the ammonia boiler is also fixed. Recall that m] is an important 
consideration in the heat transfer analysis of the steam condenser/ammonia 
boi 1 er . 
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The thermal efficiency of the ammonia cycle is: 

4 5 

iil ''Tj -jSa '*P Net Power (NH 3 ) (equation 72) 

"^th NH 3 ■ 

NH 3 NH 3 

c. The Combined Cycle 

The thermal efficiency of the combined cycle, is the total 

of the net power outputs for both cycles divided by the heat added 
to the cycle: 

^th cc = ^st ^ ^NH 3 (equation 73) 

^A steam 

recall that Pst = Qa steam x '^th ^^eam 

Pnh 3 = Qa NH 3 X nth NH 3 (equation 75) 

where NH 3 = Qa steam ‘ ’^th (equation 76) 

therefore, nth cc * ^steam “ ’isteam^^th NH 3 (equation 77) 

There are several other factors which cause the overall efficiency 
to be lower than the thermal efficiency. The first is the pumping 
requirements for the ammonia condenser discussed in the preceding 
chapter, and labelled PP. The second consideration is the power 
loss associated with converting the mechanical power of the turbines 
to electrical power, and may be expressed as Pg* the generator efficiency. 
The third factor is the power that must be expended to prepare the coal 
and to drive the auxiliary equipment in the plant. This factor is 
expressed as ria> the auxiliary efficiency. 
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The fourth and final factor is the boiler efficiency nb which is 
a measure of how much of the energy released in the combustion of 
the fuel actually is transferred to the working fluid. 

The cycle efficiency is the efficiency of the cycle including 
the condenser pumping power: 

^c = ^*^st ^NH 3 ^ ~ = ritb cc " steam (equation 78) 

% steam 

Thus the overall plant efficiency is: 

= nc X tig X rig X n |3 (equation 79) 

In this chapter, the thermodynamic analysis of the combined cycle 
was undertaken. In the next chapter, the thermodynanic and component 
considerations are combined in the economic analysis. 
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VI. ECONOMIC ANALYSIS 

The economic factors of concern in the analysis are the 

capital and operating expenses. For simplicity, only the cost 

differences between the combined cycle plant and 'Bull Run' 

will be considered. 

The significant capital expenses are: 

(!) Turbomachinery - specifically the ammonia and LP 
steam turbines. 

(2) Heat exchangers - specifically the steam condenser/ 
ammonia boiler. 

(3) Financing costs 

The significant operating expense is the fuel cost which is 
derived from the cycle overall efficiency. 

The major difference between 'Bull Run' and the combined cycle 
in turbomachinery costs is the addition of the ammonia turbine and 
the removal of the last stage or two of the low pressure steam 
turbine. 

Due to the large specific volume of steam at condensing conditions 

in conventional power plants, the last few stages of the LP steam 

1 8 

turbine are quite large and costly. Data from DeLaval indicates 
that the cost of a steam turbine varies approximately linearly with 
power output as the backpressure is increased. This data fixes the 
cost variation with power at $31 per kilowatt of mechanical power output. 
The vapor pressure of ammonia at condensing conditions is much 
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higher than that of steam. Consequently, the specific volume of 

ammonia is 100 - 300 times lower than steam, causing the ammonia 

turbine to be smaller and cheaper than the steam turbine it replaces. 

19 

The selection of $5/kw for the ammonia turbine is consistent with 
previous studies and existing NH 3 turbine costs corrected to 1976 
dollars. 

The net power output of the combined cycle under consideration 
is the same as that of the 'Bull Run' plant scaled up to 1000 MW. 
Therefore, any power lost by the steam turbine by the addition of the 
bottoming cycle must be compensated for by the ammonia turbine. 

The difference in turbine costs between the binary-cycle plant 
and 'Bull Run' may be expressed as the product of the NH 3 turbine 
power and the difference in cost per kilowatt between the steam and 
ammonia turbines. Thus, 

TURBOMACHINERY COST SAVINGS = x ($31/kw - $5/kw), (equation 80) 
where Pf^H 3 was derived as a function of the operating parameters in 
the last chapter. 

The most significant heat exchanger cost in the combined cycle is 
that of the steam condenser/ ammonia boiler. The total area of the SC/AB 
is the sum of the boiling region area, Ag,and the non -boiling area 
^NB, introduced in chapter four. 

Aj = Ag + Anb (equation 81) 

20 

Utilizing heat exchanger cost information , the cost of the SC/AB 
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per unit area is assumed to be $8/ft^. The cost of the steam condenser/ 
ammonia boiler is: 

COST SC/AB = A-j- X $8/ft^ (equation 82) 

Combining equations (80) and (82) the additional initial capital 
expense of the combined cycle is: 

INITIAL CAPITAL COST = Aj x $8/ft^ - x $26/kw (equation 83) 

In practice, the funds for capital equipment are borrowed, 
incurring finance charges. 

For a loan where the annual repayment rate is the same for every 

year, the total cost of capital equipment including finance charges is: 

TOTAL COST , INITIAL CAPITAL 1 (1 + D" ^ „ (eouation 84) 

OF EQUIPMENT COST “ + i)n. , * " (equation 



where n is the loan duration and i is the interest rate per annum. 

It must be noted that the turbomachinery cost per kw and heat 
exchanger cost per square foot of surface area represent 1976 values. 

The analysis and optimization procedure is easily modified to reflect 
price changes. 

The operating costs for the combined cycle when compared with 
those of the unmodified 'Bull Run' plant are reflected in the differences 
in fuel consumption between the two cases. The heat input to a cycle 
is the product of the fuel mass flow rate and the heating value of the 
fuel: HEAT INPUT = mj-Qg] x HHVcoal (equation 85) 

The power output is the product of the heat input and the overall 
efficiency, POWER = HEAT INPUT/noA (equation 86) 
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Combining equations (85) and ( 86 ), the coal mass flow rate is: 

= POWER OUTPUT (equation 87) 



^ coal 



nOA X ^^^'^coal 



The difference in coal flow between the 'Bull Run' plant and 
the combined cycle is: 



"'BULL RUN 



- m 



cc 



^ (POWER 0UTPUT)br . (POWER 0UTPUT)cc (equation 88 ^ 



nOA X HHVcoal IqA ^ ^^^'^coal 
BULL 
RUN 



cc 



Since the analysis has assumed that the steam portion of the 
combined cycle except for the condensing region is identical to 
'Bull Run' and that both plants have the same net power, equation ( 88 ) 
may be rewritten as follows: 

(equation 89) 



• _ POWER OUTPUT 

"’BR - "icc ■ HHV 



coal X Tia X Tig X Tib 



1 - 1 
^\)br 



cc 



where rig* hg» 15 refer to equation (79) and ti^ is the cycle 
efficiency, a modification of thermal efficiency which includes condenser 
pumping requirements. 

Equation (89) indicates that if the combined cycle efficiency > 

(n^)cc is higher than the 'Bull Run' cycle efficiency (tic)br 
be saved. 

The economic gain of reduced fuel consumption is: 

DOLLAR SAVINGS PER HOUR = COAL COST $/lbm x (liigj^ - m^^) 

Equation (90) could be used to extrapolate the savings over the plant 
lifetime, however the term 'power output' appears in equation (89). 

It would be too optimistic to design a combined cycle based on 100% 



- 65 - 



load demand over the entire plant life. The quantity load factor, 
LF, expresses the ratio of the expected average power output of the 
cycle to its maximum capacity. Equation (89) becomes: 



f" BR - ^'cc 



^XIMUM POWER OUTPUT x LF 





(eq. 91) 



Then the fuel savings in dollars over the plant lifetime is: 



LIFETIME 

FUEL SAVINGS COAL COST $/Ton x ’' 200 ^ 



(eq. 92) 

X 1 Ton X 8760 hr x 
W 



where n is the expected plant lifetime in years, as well as the loan 
duration. 

Combining the operating and capital expenses into one equation, 
the life cycle benefit in dollars of the combined cycle over the 
unmodified 'Bull Run' plant is: 

SAVINGS FUEL SAVINGS OF EQUIPMENT tguMiiuiN 

The purpose of the optimization procedure is to maximize the 
'NET SAVINGS' of equation (93), and determine the operating parameters, 
in terms of the three independent variables T 5 , TF(1) and TF(5), which 
yield the maximum economic benefit. 



n yrs 
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VII. OPTIMIZATION PROCEDURE 

In this chapter, the procedure for the selection of the 
optimum design bottoming cycle based on economic and engineering 
factors will be outlined. 

The optimum design may be specified in terms of three operating 
parameters previously identified as the independent variables in the 
analysis, they are: 1 

(1) The steam condensing temperature, T 5 . 

(2) The ammonia boiling temperature, TF(1). 

(3) The ammonia condensing temperature, TF(5). 

Once values for these variables are chosen, a thermodynamic 
analysis of the combined cycle is undertaken which yields the combined 
cycle thermal efficiency and the power split as the output. Power 
split is defined as the amount of power produced by each of the two 
cycles. ’ 

Next, the heat transfer analysis for the steam condenser/ammonia?' 
boiler specifies the size of the heat exchanger. The ammonia condenser 
is then analyzed with the condenser pumping power as the result. 

The information derived from these analyses is combined in 
equation (93), which combines capital and operating costs into a value of 
net savings over the plant lifetime. 

A flow diagram of the optimization procedure is shown as figure 
17. Note that the value of net savings resulting from the analysis is 
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FIGURE 17 OPTIMIZATION PROCEDURE FLOW CHART 
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very sensitive to the economic factors used in the cost computations. 

The most significant factors are fuel and heat exchanger costs. 

Thus far, the optimization has produced a value of net savings 
of the combined cycle plant over that of the 'Bull Run' plant scaled 
up to 1000 MW for one set of the three variables T^, TF(1), and TF(5). 

In order to determine which set will yield the optimum plant, the entire 
analysis must be carried out for all possible combinations of the three 
independent parameters within the range of values where the bottoming 
cycle is expected to operate. 

The analysis may be accomplished independent of environmental 
conditions with regard to selection of values of the steam condensing 
temperature, T$, and the ammonia boiling temperature TF( 1 ). However, 
the selection of the ammonia condensing temperature, TF(5) is very much 
affected by the temperature of the condenser coolant water at the inlet. 

As discussed in chapter four, the relationship between the condenser 
water inlet temperature and the ammonia condensing temperature is restricted by 
considerations of pumping power and tube corrosion resulting from high 
water mass flow rates and velocities. 

Utilizing equations (13) and (78) the cycle efficiency may be 
calculated for various values of temperature difference between the con- 
densing ammonia and the inlet water for fixed parameters T 5 , TF(1), 
and TF(5). Additionally, the water velocity through the tubes as a func- 
tion of the same parameters may be calculated. The results of these 
calculations are plotted as figure 18. The values of steam condensing 
temperature, ammonia boiling temperatures, and ammonia condensing tem- 
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VELOCITY 
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FIGURE 18 CONDENSER WATER VELOCITY AND CYCLE EFFICIENCY VS. DIFFERENCE 
BETWEEN NH 3 CONDENSING TEMPERATURE AND WATER INLET TEMPERATURE 
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perature used for figure 18 are representative, near optimum values. 

Notice that the cycle efficiency is at a maximum at approximately a 

15°F temperature difference between the condensing ammonia and the 

inlet coolant water. Furthermore, the same temperature difference 

yields a water velocity of 12.5 ft/sec, close to the corrosion limit 

of 15.0 ft/sec. Thus it appears that the optimum selection of TF(5), 

the ammonia condensing temperature is 15°F above the available water 

source temperature derived from environmental information. This conclu- 

21 

sion agrees with earlier condenser design studies. 

Replacing the ammonia condensing temperature TF(5) as one of the 
independent parameters by the water inlet temperature Tjj^ and keeping 
the 15°F difference between TF(5) and Tji^ for all possible designs, the 
optimization computations for different sets of the independent par- 
ameters may now proceed. 

Since the combined cycle under consideration is compared to the 

'Bull Run' plant, the environmental conditions of the 'Bull Run' plant 

are used. At the 'Bull Run' site, the average condenser inlet temper- 
22 

ature is 55°F , and this value is incorporated into the optimum design 

analysis. However, other values of condenser inlet temperature are 
included in the results. 

Therefore, only the steam condensing and ammonia boiling tem- 
peratures need be varied to find the optimum design. In practice, 
the steam condensing temperature is fixed and the temperature difference 
between the condensing steam and boiling ammonia is varied. To present 
the results on a non-dimensional basis, this temperature difference is 
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expressed as the difference between the steam condensing temperature 
and the entropy averaged temperature of heat addition to the ammonia, 

T, divided by the steam condensing temperature, 

AT/T^ = (T^ - T)/ Tj (equation 94) 

where T is a function of the ammonia boiler inlet and outlet conditions 

In order to carry out the optimization procedure, the many 
engineering and economic factors previously discussed must be given 
values. The optimization computer program (see appendix one), will 
carry out the analysis for any values of these factors. This flexi- 
bility enables this optimization to be used for future bottoming cycle 
design utilizing up-to-date economic considerations. The symbols and 
values used in this analysis are shown in table 4. 

Thus far, the derivation of the analysis and the procedure for 
identifying the optimum design combined cycle have been presented. 

In the next chapter, the results of the optimization calculations are 
presented along with the specifications of the optimum bottoming cycle 
compatable with the 'Bull Run' plant. 
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EFFICIENCIES 


TEXT 

SYMBOL 


COMPUTER 

SYMBOL 


VALUE 


1. 


Auxiliary Efficiency 


’^a 


NA 


94 % 


2. 


Boiler Efficiency 


% 


R4 


89 % 


3. 


Bull Run Cycle Efficiency 




R5 


45.78% 


4. 


Generator Efficiency 


"g 


R7 


99 % 


5. 


Pump Efficiency 


"P 


NP 


90 % 


6. 


Turbine Isentropic Efficiency 


>^is 


NS 


84 % 


7. 


Turbine Mechanical Efficiency 
(Includes exit, stage and 
mechanical losses) 


^M 


NM 


97 % 




HEAT EXCHANGER SPECIFICATIONS 








8. 


Ammonia Boiler Tube Diameter 


D 


D2 


1" f 


9. 


Ammonia Boiler Tubes in Vertical 
Bank 


n 


NOT 


70 


10. 


Ammonia Condenser Tube Diameter 


D 


D3 


.875" 


n. 


Ammonia Condenser Tubes in Vertical 
Bank 


n 


INCLUDED IN 
A CONSTANT 


200 




ECONOMIC FACTORS 








12. 


Coal Cost (111. #6) 


- 


R6 


$50/Ton 


13. 


Heat Exchanger Cost 


- 


R1 


$ 8/ft^ 


14. 


Higher Heating Value (111. #6) 


HHV 


RF 


10788. BTU/lbm 


15. 


Interest Rate 


i 


IR 


12% 


16. 


Load Factor 


LF 


R8 


.8 


17. 


Loan Duration and Plant Life 


n 


R9 


30 years 


18. 


Turbine Cost (Ammonia) 


« 


R3 


$ 5/kw 


19. 


Turbine Cost (LP Steam) 


• 


R2 


$31/kw 



TABLE 4 SUMMARY OF ENGINEERING AND ECONOMIC FACTORS 
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VIII. OPTIMIZATION RESULTS 

Utilizing the optimization procedure outlined in the last 
chapter, the thermal efficiency and lifetime net savings for each 
set of combined cycle parameters are calculated. Assuming the 
average condenser water inlet temperature to be 55°F, and 
incorporating the engineering and economic factors enumerated in 
table 4, an optimum design may be predicted. 

Figure 19 is a plot of the combined cycle efficiency versus 
the non-dimensional temperature difference AT/T 5 for several different 
values of Ts, the steam condensing temperature. As expected, the cycle 
efficiency increases as the temperature difference decreases. 

Depending on the value of AT/T 5 , either a steam condensing temperature 
of 131°F or 147°F yields the maximum efficiency. This efficiency peak 
is caused by wetness considerations discussed in chapter 4. 

However, as the non-dimensional temperature difference decreases, 
the area of the steam condenser/ammonia boiler increases. The 
opposing trends of efficiency and heat exchanger size variation with 
the non-dimensional temperature difference are the most significant factors 
in the optimization routine. Figure 20 shows the steam condenser/ 
ammonia boiler areas versus the non-dimensional temperature difference 
for several values of T^, the steam condensing temperature. Note that 
the heat exchanger area is slightly smaller for T^ = 147°F than for 
131°F. This trend is caused by the change of the steam properties which 
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NON-DIMENSIONAL TEMPERATURE DIFFERENCE, AT / Tg 



FIGURE 19 COMBINED CYCLE EFFICIENCY VERSUS NON-DIMENSIONAL 
TEMPERATURE DIFFERENCE FOR VARIOUS VALUES OF 
STEAM CONDENSING TEMPERATURE. WATER INLET 
TEMPERATURE = 55°F. 



- 75 - 




NON-DIMENSIONAL TEMPERATURE DIFFERENCE, aT / T 5 



FIGURE 20 STEAM CONDENSER/ AMMONIA BOILER HEAT TRANSFER AREA 
VERSUS NON-DIMENSIONAL TEMPERATURE DIFFERENCE FOR 
VARIOUS VALUES OF STEAM CONDENSING TEMPERATURE. 
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affect the steam side heat transfer coefficient (seeequation (23)) as 
the condensing temperature varies. 

Combining the effects of cycle efficiency, heat exchanger 
size, and the other economic factors in the cost equation (equation (93)), 
a plot of lifetime cycle savings is presented as figure 21. 

It is obvious that the maximum savings occurs where the steam 
condensing temperature = 147°F, and the non-dimensional temperature 
difference AT / T5 = .09. This value of AT / T5 corresponds to an 
ammonia boiling temperature TF(1) = 142°F. 

The predicted lifetime combined cycle savings over the 'Bull Run' 
plant is $98 million or about .047 cents/kw - hr. The combined cycle 
efficiency is 47.2%, about 1.4% higher than the 'Bull Run' plant 
with an initial additional capital cost of only $9.1 million. 

The flow chart and thermodynamic specifications of the combined 
cycle are presented as figure 22 and table 5 respectively. 

In order to determine the combined cycle performance in the off- 
design condition, several assumptions must be made. They are: 

(1) Maintain steam condensing and ammonia boiling 
temperature at design points. 

(2) Maintain a 15°F temperature difference between the 
anmonia condensing temperature and the coolant water 
inlet temperature. 

(3) Utilize DeLaval estimate of off-design ammonia turbine 
performance (figure 23). 

Carrying out the off-design analysis, incorporating the exit 
loss criteria from figure 23, a plot of combined cycle efficiency versus 
coolant water inlet temperature is shown as figure 24. For water inlet 
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FIGURE 21 PLANT LIFETIME NET SAVINGS VERSUS NON-DIMENSIONAL 
TEMPERATURE DIFFERENCE FOR VARIOUS VALUES OF 
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TEMPERATURE = 55°F 
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FIGURE 22 COMPONENT FLOW CHART FOR COMBINED CYCLE 




a 



STREAM 


NO. 


DESCRIPTION 


P(psi) 


T(°F) 


h(BTU/lbm) 


m( 10 ^ Ibm/hr) 


Fuel 


a 


Raw fuel (13% H 2 O) 


14.7 


60 


HHV 10788 


.81 




Fuel 


a' 


Prepared fuel 
( 6 % H 2 O) 


14.7 


60 


HHV 10788 


.75 




Air 


b 


Combustor inlet air 


14.7 


60 


124.3 


8.21 




Exhaust c 


Products of 
combustion 


14.7 


300 


189. !• 


8.82 




Ash 


d 


Ash 


14.7 


- 


- 


.08 




H 2 O 


e 


Water extracted 
from fuel 


14.7 


60 


28.0 


.06 




H 2 O 


1 


Condenser exit 


3.45 


147 


115.0 


4.68 




H 2 O 


2 * 


Boiler inlet 


4210 


550.1 


544.9 


7.31 




H 2 O 


3 


Boiler exit/HP 
turbine inlet 


3515 


1000 


1420.8 


7.31 




0 

CVJ 


4 


HP turbine exit/ 
reheater inlet 


600 


551.7 


1255.4 


5.58 




H 2 O 


5 


Reheater exit/IP 
turbine inlet 


540 


1000 


1519.6 


5.17 




H 2 O 


6 


IP turbine exit/ 
condenser inlet 


3.45 


147 


1078.7 


4.68 X = 


.954 


NH 3 


7 


NH 3 condenser exit 


128.8 


70 


121.1 


8.35 




NH 3 


8 * 


NH 3 boiler inlet 


390.3 


124 


182.4 


10.0 




NH 3 


9 


NH 3 boiler exit 


390.3 


142 


632.4 


10.0 




NH 3 


10 


NH 3 turbine exit 


128.8 


70 


612.2 


8.35 X = 


.967 


H 2 O 


11 


Coolant water inlet 


14.7 


55 


23.0 


513 




H 2 O 


12 Coolant water exit 
* INCLUDES FEEDHEATING 


14.7 


63 


31 .0 


513 




TABLE 


5(PART 1) THERMODYNAMIC SPECIFICATIONS OF 


THE COMBINED 


CYCLE 





(PART 2 ON NEXT PAGE) 



- 80 - 



COMPONENT POWER 

Steam boiler 
Steam reheater 

HP steam turbine 272.2 MW 

IP steam turbine 682.6 MW 

Steam condenser/NH 3 boiler 
Ammonia turbine 123.9 MW 

Ammonia condenser 



Gross total power 
Generator losses 
Auxiliary power 
Condenser pumps 
NET POWER 



1078.7 MW 
- 10.8 
-64.1 
- 3.8 
1000.0 MW 



HEAT TRANSFER 

1876.0 MW 
400.5 MW 

1321.7 MM 
1201.9 MW 



Combined cycle: Thermal efficiency = 47.38% 

Cycle efficiency = 47.2% 

Overall efficiency = 39.1% 

Steam cycle : nth " 41.94% 

Ammonia cycle : 9-37% 



TABLE 5(PART 2) THERMODYNAMIC SPECIFICATIONS OF THE COMBINED CYCLE 
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FIGURE 23 EXIT LOSSES AMMONIA TURBINE IN OFF DESIGN CONDITIONS 



49.0 




Tin ('’f) 

FIGURE 24 OFF-DESIGN PERFORMANCE OF COMBINED CYCLE 
-83- 






I 



temperature above 55° F, the pressure ratio used in figure 23 becomes 
greater than one. For pressure ratios greater than unity, the exit 
loss is assumed to be 1%. 

Optimization results for condenser water inlet temperatures of 
35°F, 45°F, and 65°F are shown in figure 25 through 30. The parameters 
of the optimum design for each value of condenser inlet temperature are 
shown in table 6. 



CONDENSER 

INLET 


STEAM 

CONDENSING 


NH3 

CONDENSING 


NHo 

BOILING 


NET 


TEMPERATURE 


TEMPERATURE 


TEMPERATURE 


TEMPERATURE 


SAVINGS 


35°F 


131°F 


50°F 


126°F 


$233 million 


45°F 


147°F 


60°F 


142°F 


$162 million 


55°F 


147°F 


70°F 


142°F 


$ 98 million 


65°F 


147°F 


80° F 


142°F 


$ 29 million 



TABLE 6 OPTUMUM DESIGN PARAMETERS 



The net savings for each optimum design taken from table 6 
is plotted in figure 31. As expected, the combined cycle is most 
attractive where the ambient temperature is lowest. More significantly, 
figure 31 indicates at what value of condenser inlet temperature the 
addition of a bottoming cycle is no longer economically viable, approximately 
69°F. 

In this chapter, information has been presented to permit the 
selection of an optimum ammonia bottoming cycle with the available water 
source temperature as the only input variable. 
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NON-DIMENSIONAL TEMPERATURE DIFFERENCE aT / Ts 

FIGURE 25 COMBINED CYCLE EFFICIENCY VERSUS NON-DIMENSIONAL 
TEMPERATURE DIFFERENCE 
CONDENSER INLET TEMPERATURE = 35°F 
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250 



Tin = 35°F 
TF(5) = 50°F 



200 .. 



LIFETIME NET 
SAVINGS ($ X lOO 



T$ = 131 




150 



NON-DIMENSIONAL *¥JmPERATURE DIFFERENCE / Is 

FIGURE 26 COMBINED CYCLE NET SAVINGS VERSUS NON-DIMENSIONAL 
TEMPERATURE DIFFERENCE 
CONDENSER INLET TEMPERATURE = 35°F 
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NON-DIMENSIONAL TEMPERATURE DIFFERENCE AT / Ts 

FIGURE 27 COMBINED CYCLE EFFICIENCY VERSUS NON-DIMENSIONAL 
TEMPERATURE DIFFERENCE 
CONDENSER INLET TEMPERATURE = AB'^F 



200 



T 



IN 



45°F 



TF(5) = 60°F 




50 



+ 



+ 



.01 



.02 



NON-DIMENSIONAL TEMPERATURE DIFFERENCE AT/Ts 



FIGURE 28 COMBINED CYCLE NET SAVINGS VERSUS NON-DIMENSIONAL 
TEMPERATURE DIFFERENCE 
CONDENSER INLET TEMPERATURE = 45°F 
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NON-DIMENSIONAL TEMPERATURE DIFFERENCE AT / Tj 

FIGURE 29 COMBINED CYCLE EFFICIENCY VERSUS NON-DIMENSIONAL 
TEMPERATURE DIFFERENCE 
CONDENSER INLET TEMPERATURE = 65°F 
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= 65°F 




FIGURE 30 COMBINED CYCLE NET SAVINGS VERSUS NON-DIMENSIONAL 
TEMPERATURE DIFFERENCE 
CONDENSER INLET TEMPERATURE 65°F 

- 90 - 



250 



200 



150 



LIFETIME NET 
SAVINGS 

($ X 10°) 100 



50 



0 




35 45 55 65 75 

CONDENSER INLET TEMPERATURE (°F) 

FIGURE 31 OPTIMUM DESIGN NET SAVINGS VERSUS CONDENSER 
INLET TEMPERATURE 



- 91 - 



IX. SUMMARY AND CONCLUSION 

a. SUMMARY 

The purpose of a bottoming cycle is to improve the efficiency 
of a power plant by enabling the temperature of heat rejection to be 
lower than what is practical for a steam plant. 

By improving the cycle efficiency, a bottoming cycle equipped 
power plant will consume less fuel, reject less heat to the environment, 
and deposit fewer harmful products of combustion into the air. 

Ammonia seems to be the best choice as the working fluid for a 
bottoming cycle designed to operate in sub-position to a conventional 
steam plant such as 'Bull Run'. 

The bottoming cycle is defined by three operating parameters: 

(1) T 5 , the steam condensing temperature. 

(2) TF(1), the ammonia boiling temperature. 

(3) TF(5), the ammonia condensing temperature. 

In order to determine the operating parameters of the optimum 
bottoming cycle, a thermodynamic analysis was conducted to find the cycle 
efficiency as a function of the operating parameters. Next, the steam 
condenser/ammonia boiler and turbomachinery sizes were calculated. 

These operating and capital expenses are combined into an equation which 
yields the net lifetime savings of a combined cycle plant over the 
unmodified steam plant. The set of operating parameters which gives the 
maximum economic benefit is taken as the optimum design. 
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The steam condensing temperature which gives maximum efficiency 
is that value which minimizes the total power lost due to wetness in 
both the steam and ammonia cycles. 

Limitations on the ammonia condensing temperature are due to cor- 
rosion and pumping power considerations on the water side of the condenser. 
The ammonia condensing temperature for optimum performance is 15°F 
above the condenser water inlet temperature. 

Thus, the economic analysis hinges on the selection of the ammonia 
boiling temperature, or more precisely the difference between the 
steam condensing temperature and the entropy averaged temperature of 
heat addition to the ammonia. As this temperature difference decreases, 
the thermodynamic efficiency of the combined cycle increases but the 
area of the steam condenser/ammonia boiler increases. The economic 
analysis indicates what value of this temperature difference makes the 
addition of an ammonia bottoming cycle most attractive. 

Given the present economic factors enumerated in table 4, the 
optimum bottoming cycle to be installed at the 'Bull Run' site has 
the following operating and economic parameters: 



(1) 


Steam condensing temperature 


Ts = 147°F 


(2) 


Ammonia boiling temperature 


TF(1) = 142°F 


(3) 


Ammonia condensing temperature 


TF(5) = 70°F 


(4) 


Condenser water inlet temperature 


Tin = 55°F 


(5) 


Cycle net power 


P = 1000 MW 


(6) 


Overall cycle efficiency 


riQA ” 39.1% 
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I 




(7) Cycle efficiency increase over An = 1.4% 

•Bull Run' 

(8) Lifetime net savings $98 million or .047^/kw-hr 

(9) Added initial capital cost $9.1 million 

b. CONCLUSION 

The ammonia bottoming cycle is economically and technologically 
feasible using existing technology and based on current economic 
conditions, provided the condenser water inlet temperature is less 
than 69°F. If the cost of fuel continues to increase, which is likely, 
the bottoming cycle becomes even more attractive. 

The addition of a bottoming cycle seems to be suitable for large 
conventional utility steam plants, and pressuri zed-water type nuclear 
plants employing water-cooled condensers. However, due to the large 
size of the steam condenser/ammonia boiler which is on the order of 
1.5 X 10® square feet of surface area, the bottoming cycle is not deemed 
practical for applications where power plant size is a major consideration. 
Additionally, the possibility of a major ammonia leak precludes the 
installation in a closed-in environment such as a ship. 
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APPENDIX I THE COMPUTER PROGRAM 



a. INTRODUCTION - The computer program used in the optimization/ 
analysis was run on the M70 and M80 computers at the joint 
Mechanical/Civil Engineering computer facility at MIT. The programming 



language used was Fortran IV. 



b. 

VARIABLE 

A 

AB 

ANB 

AT 

A2 

A5 

A 6 

A7 

BETA 

B1 

C 

CAPT 

CP 

Cl 

C3 

C4 

C5 



IDENTIFICATION OF COMPUTER VARIABLES 
REPRESENTATION 

Intermediate variable which sums AB's 

Area of one quality region of boiling section in SC/AB 

Area of non-boiling section of SC/AB 

Total area of steam condenser/NH 3 boiler 

NH 3 condenser area 

Difference between density of saturated liquid and 
saturated vapor of NH 3 at condensing temperature 

Density of saturated liquid NH3 at condensing temperature 

Difference between NH3 condensing temperature and the 
tube wall temperature 

Part of condensing heat transfer coefficient 

Average quality in a section of the boiling region of 
NH 3 boiler 

Heat transfer coefficient on boiling side of SC/AB 
(first approximation) 

Additional capital cost of bottoming cycle 
Specific heat at constant pressure of NH 3 

Heat transfer coefficient on boiling side of SC/AB( final value) 

Heat transfer coefficient on NH 3 side of non-boiling section 
of SC/AB 

Heat transfer coefficient on steam side of non-boiling 
section of SC/AB (first approximation) 

Error between C3 and C4 in iteration loop 
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VARIABLE 

C6 

D 

DELTA 

DP 

D1 

D2 

D3 

E 

EFF 

EFFOA 

EFFST 

ENTH 

ETA 

El 

E2 

E3 

F 

FI( ) 
FII( ) 

G 
GU 
H( ) 

HC 

HCOND 
HF( ) 
HFG( ) 



REPRESENTATION 

Final value of steam condensing heat transfer coefficient 

Heat transfer coefficient on condensing side of SC/AB 
(first approximation) 

Difference between the steam condensing temperation (T$) 
and the average temperature of NH3 heat addition (TBAR) 

Pressure drop on water side of NH3 condenser 

Heat transfer coefficient on condensing side of SC/AB 
(final value) 

NH3 boiler tube diameter 

NH3 condenser tube diameter 

Error in heat flux calculation in SC/AB (C-D) 

Thermal efficiency of NH3 cycle 

Thermal efficiency of combined cycle excluding condensing 
pump requirements 

Efficiency of steam cycle at specified steam condensing 
temperature 

Corrected enthalpy used in NH3 condensing calculation 
Overall thermal efficiency of combined cycle 

Absolute difference between first and second approximation 
in NH3 condenser water outlet temperature routine 

Error in first approximation of heat flux in NH3 condenser 
Absolute value of E2 
Chen factor for boiling 

Intermediate function used to calculate HG( ) 

Intermediate function used to calculate SG( ) 

NH3 mass flow per unit area in NH3 boiler 

Mass flow rate/unit area of water in NH3 condenser 

Enthalpy (NH3) 

Convective component of the heat transfer coefficient in 
the NH3 boiler 

Condensing NH3 heat transfer coefficient 
Enthalpy of saturated liquid (NH3) 

Latent heat of evaporation (NH3) 
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VARIABLE 
HG( ) 
HNUB 

HS( ) 

HW 

I 

IR 

KF 

KFl 

KW 

L 

LMTD 
ME( ) 
MNH 3 
MW 
N 

NM 

NOT 
NP 
NS 
P( ) 

PF 
PG 
PH20 
PNH3 
PP 
PUMPW 
PX( ) 

Q 

QA 



REPRESENTATION 

Latent heat of evaporation vapor (NH^) 

Nucleate boiling component of heat transfer coefficient 
in NHg boiler 

Enthalpy (based on isentropic expansion in NH3 turbine) 

Convective heat transfer coefficient of water side 
of NH^ condenser 

Counting integer in DO LOOP 

Capitalization rate 

Thermal conductivity of NH^ in NH3 boiler 
Thermal conductivity of NH^ in NH^ condenser 
Thermal conductivity of water 
Temperature step side in AT variation routine 
Log mean temperature difference SC/AB 
NH^ turbine mass flow extraction for feedheating 
Mass flow rate through NH3 boiler 
Mass flow rate through NH^ condenser 
Counting integer in DO LOOP 

Turbine efficiency after mechanical, exit and stage 
losses 

Number of tubes in NH^ boiler 
Pump efficiency 

Isentropic efficiency of turbines 

Pressure (in thermodynamic properties of NH3 sub routine) 

Density of NH3 saturated liquid 

Density of NH^ saturated vapor 

Power produced by steam cycle 

Power produced by NH^ 

NH3 condenser water pumping power 

Total NH^ feed pumping power 

Intermediate function used to calculate P( ) 

Heat added to NH cycle 
3 

Enthalpy added to NH^ in SC/AB 



- 101 - 



VARIABLE 

QAD 

QB 

QNB 

QR 
RETP 
R1 
R2 
R3 
R4 
R5 
R6 
R7 
R8 
R9 
S( ) 

SAVE 

SF{ ) 

SG( ) 

ST 
SI 
T{ ) 

TAV 
TBAR 
TF( ) 

TIN 

TOUT 

TOUTB 

TS 

TURBW 



REPRESENTATION 

Heat added to steam cycle 

Total heat added to boiling section of NH^ boiler 

Total heat added to non-boiling section of NH^ boiler 

Total heat rejected by combined cycle 

Two-phase flow Reynolds number 

2 

Heat exchanger cost $/FT 
Steam turbine cost $/MW 
NHj turbine cost $/MW 
Boiler efficiency 
Bull Run thermal efficiency 
Coal cost $/T0N 
Generator efficiency 
Load factor 

Plant life and loan duration 
Entropy 

Lifetime savings gained by bottoming cycle equipped 
plant over base case 

Entropy of saturated liquid (NH3) 

Entropy of saturated vapor (NH^) 

Surface tension of NH3 in NH3 boiler 

Boiling suppression factor in two-phase flow 

Temperature (*‘R) 

Average water temperature in NH^ condenser 
Average temperature of heat addition of NH3 
Temperature (°F) 

Inlet water temperature in NH^ condenser 

Outlet water temperature in NHo condenser (first approx- 
imation ) 

Outlet water temperature in NH^ condenser (second approx- 
imation) ^ 

Steam condensing temperature (°F) 

Total NH^ turbine power 
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■m 




VARIABLE 


REPRESENTATION 


TT 


NH^ condenser tube wall temperature 


TW 


SC/AB tube wall temperature 


U 


Overall heat transfer coefficient 


UB 


Overall average heat transfer coefficient in boiling 
section of the SC/AB 


UC 

UF 


Overall heat transfer coefficient in NH^ condenser 
Absolute viscosity of NH3 saturated liquid in NH^ boiler 


UFl 


Absolute viscosity of NH3 saturated liquid in NH^ condenser 


UG 


Absolute viscosity of NH^ saturated vapor in NH^ boiler 


UNB 


Overall heat transfer coefficient in non-boiling section 
of the SC/AB 


UW 


Average absolute viscosity of water in NH3 condenser 


VF( ) 
VFG( ) 
VG( ) 
VW 


Specific volume of NH saturated liquid 
VG( ) - VF( ) ^ 

Specific volume of NH3 saturated vapor 
Velocity of water in NH3 condenser 


WP( ) 


NH3 feed stage pump work 


WT( ) 
X( ) 
xs( ) 

XTT 


NH3 turbine stage work 
Quality 

Quality based on isentropic expansion in NH3 turbine 
Martinelli parameter for boiling 


YTT 


Inverse of XTT 


Y 


% quality in NH3 boiler section (B1 x 100) 


Z 


Delta/Ts(°R) 
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c. PROGRAM DESCRIPTION 



For clarity, the computer program on the following pages is 
divided into nine zones, labelled 'A' through 'I'. Each zone 
represents a distinct function in the optimization routine. 

Zone 'A' includes the identification of variables as either 
integers or real numbers. The 'dimension' statements reserve core 
space for the subscripted variables. 

Zone 'B' consists of the engineering and economic factors 
listed in table four. It is these factors which the economic 
analysis hinges upon. Accordingly, these parameters are easily 
changed to allow an up-to-date analysis. Also included in zone 
'B' are the ammonia condensing temperature and the coolant water 
inlet temperature for the plant under consideration. 

The steam condensing temperature, steam cycle efficiency, 
and part of the steam side heat transfer coefficient are included 
in zone 'C. Also included are the temperature difference between 
the condensing steam and boiling ammonia, as well as the pump and 
turbine efficiencies. 

Part 'D' establishes the three intermediate temperature levels 
in the ammonia cycle and generates the thermodynamic properties of 
ammonia needed for the analysis. The equation of state for the ammonia 
are empirical curves, fitted to the tabulated data^^- 
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The thermodynamic analysis of the ammonia cycle and the 
calculation of the combined cycle thermal efficiency occurs 
in zone 'E'. It is in this section where the non-dimensional 
temperature difference between the condensing steam and boiling 
ammonia is calculated. 

In zone 'F' the heat transfer properties of the ammonia are 

25 

calculated from empirical relations. 

Zone 'G' is the most complex part of the entire program. In 
this section the steam condenser/ ammonia boiler area is calculated. 

Note that in statements 55 and 77 the tube wall temperature is 
incrementally changed. As stated in chapter four, the heat flux 
equation for both sides of the steam condenser/ ammonia boiler can only 
be solved by interating until the wall temperature which satisfies 
the equality is found. In order to minimize the running time of the 
program, an exact equality of both sides of equation (34) is not 
demanded, but rather it is required that the two results are within 
1% of each other. Without this relaxation of the solution, the 
temperature increment required to obtain a converging iterative solution 
would increase the running time and cost of the program by several orders 
of magnitude. Note also that the heat exchanger analysis is repeated 
for the 10 quality regions discussed in chapter four, as well as for the 
non-boiling region. 

The ammonia heat transfer properties and condenser pumping power 
are computed in zone 'H'. Once again the tube wall temperature is found 
by iteration. 



- 105 - 



In zone 'I', the condenser pumping power is incorporated into 
the cycle efficiency. The initial capital, operating, and lifetime 
costs are also calculated in this region. The output statements com- 
prise the final section of the program. 
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1b HI( «) = p(r. )*(9b2.-2.964*p( v ) + .0C62E6*P( « ) * *2-4 . f9E-C6*P( N ) ** 3 ) * 



1 (-1. '^-C6 ) 

SG(N )s.6 l9b46*AlOC10(T(N ))-(?( V))/T(N) *( (1.768 3F 07/T(N)**3) + 
1 (7.7474 ?7*2.3421E 2 6* P( N ) ) /T( .< ) * * 1 1 +( 6 . 908 E 37*F(N)**5)/ 
2T(b)**19 )-. 20877 23*A L0Q1 0(P( 11 )) + ( 3. 16094 F-04*T(N))-3 3,048/ 
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lit ■ ) <-.0/ 1 + r II ( . ) 

VI ( *' ) = ( . r,' !;< o^'.a ♦ . 70 7 3 i - 0 3 • r/;'HT ( 27 1 . 4-Tf' ( N ) ) -7 . 73 '■ t-o S * 

1(271 .4-7' (.)))/< 1. + .31663*n', H(271.4-TF<1i)) + 8.fS44--03* 

. (.^71 .4-7 x( ) ) ) 

Vt'•,(^■)=Vf;(^)-V^ {!■•) 

.ii ( ;)=v,. 7*1.1 )•>•( iK(;:)-30.) SECTION D 

3*( » i3('i V(I )-17C. )-KTF( ‘0-17r. ))*.?• 

hr..( •; )=HL.< )-H' ( *. ) 

IM N ) = ( 4 )-(Hr':( M/T( N) ) 

94 C0r;iI^'Ji’ 

N - /' 

C IH- IS 

rc> 111 ,1 

)=Vr (N )* n ( 1)-r ( S))/4 .*( 14 4./77P . )*( 1 ./VP) 

-( 1 ) = :.;( 1 ) 

.M( i)=}v;( 1) 

<(1)=1. 

xs ( ) = ( 3 ( V - 1 )-y:- ( ^ ) ) /( SG ( N ) -s^( V ) ) 

HS(V )^WF(^ )4XS(:.)*:UT( 0 

j:( )-(p(N-l ).4S(V ))*‘S 

>:c?) = (r:( N)-}iF(i. ))/.‘iKC(;o 

1 ('<)=SK( N)^X(V)MSu( N)-SF(N)) SECTION E 

111 Cr.NTIM.-I 

M ) = Cif(;)-Hf(3)-WP(^)) /(f?(2)-HF(3)-kPO)) 
r )=(r:M3)-(1.-ii ( 1) )♦ (Hr (4)>VP(4) ) )/( H( 3 ) -Hf ( 4 ) -WP ( 4 ) ) 

: K( i)=(i t (4)-( 1.-yK( 1)-rL(?))*(PF(^.)^W?(^.) ))/(H( 4)-Hf(S)-Wf (5)) 

M 1)-4i (i;)-.yp( V) 

pn = (-/)4(1 ,-1^ / ( 1 ))*WP( 3)^<1.-''^( D-.'^FC 2) )*WI (4 

1 ( 1 . - 'f ( 1 )-":K(2)-r^(3 ))*WM i ) 

:: = i 

bu 200 N- 1,4^1 

Xl(^ ) = (P(N)-H( ^4 1 ))*NI- 

20c COUTIMiF 

T0/^k» = ’(.T(1)4(1.-NM1 ))*WT(2)4(1.-f'r (1)-^F(2))*WT(3)4 
1 (1 1 )-vv. c^).mK(3 ))*XT(4) 

KFF=(TlR».X-PU^Pik‘)/;i 

1. V-< = (T> (1)*PFG(1)4( (TF( 1)47:'(2))/7.)*<HF(1)-(MF(2)4^?(2))))/ 
1(ii:.( 1)-H^(;)-W^(2)) 
l>«-* i.TA='1S-T- AP 

r.r f C A= i- 1 F ST4 ( 1 . -r Fl S 1 ) FF 

L’^ :’b=( .2S*(TF( 1)-TF( 5 ) ) ) /A I CG ( ( TS- . 75*TF ( 1 ) - . 2F* TF ( S ) ) / 
KTS-IFd))) 

2 = I;-LTA/(T344S0.7) 

T 

vA')=(U'0C.*J.413I Ot )/tfFOA 
', = vAl'*< 1 .-FFF5.T) 

f'N"3 = C/(^(1)-HF( 2 2)) SECTION F 

K=('» .•I'NH 0/(1 .141 »x9‘90C0. ) 

UK 1=.2<-6-.00 137K,*<n (6)-107. ) 

Ff 1=.2M-.0C;C74H*(7f (S)-97.) 

l.T--^12.11 -.079 (*(TF( 1 

'Ji ’1.1 (.‘j*.(j0 19 3*(TF( D-9 7. ) 

U .= . J2K7 +6 .484i--C0*( IF ( 1 )-97. ) 
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•Jr=.>i'f-.001 3?“;*(T''( D-107. ) 

(u=1./VG(1) 

1V = 1./VK1) SECTION F 

Ki = .2t3-.00074o*( rr( D-77.) 

Z NiM PCILi^,.- CO.'T. UTAXlONS 

.fc866»( Kt'»» .79 )♦ (Cl •* .45)* ( t'f **. U9 )• ( ( ‘ 1 C( 1 ) )•• ,51 ) )/ 
1((:/r**.5)*(0F**.:9)Mpr,*».;u)*((VFC(1))**.75)»((T( 1))**.75>) 
HC = .02 3* (Ft'/t2)*( C;*D7/UF)**.3)*((CP*UF/KF)»‘.4) 

TW=crK( 1 )+TK(i))/2. 

55 Trt^Trt^.OI 

I'c (TW.o7.TS) Cn 70 3C7 
C3-HC*(TW-(TI ( 1)+TF( 2) )/7. ) 

C4 = FFTA* (TS-T»()**.75 
C5 = A2.3(C3-t4) 

I.- (CS.CT.15.) GO TO 55 
C6=C4/(ri-T-.’) 

;=(hc»c6)/(nc+ce) 

.4-=,,* (Hi ( 1)-HF( 2)-WH2) )/(!'( 1)-HF(2 )-KF(2)) 

AN '=v-S?/(UNn*(TS-(TF(1) + TF(2) )/2. )) 

A = 0 , 

1=5 

nc tfH 1 = 5, 95, 1 C SECT I ON 6 

Y = F1.CAT( I ) 

M = Y/10C . 
b=Y/5. 

X7T=(( ( 1 1 )/: 1)**. 9)*( (PC/PF)**.* )*( CJF/UG)**. 1) 

YT1=1./XTT 

IF (YT7.CT.1.) of) TO *6 
f = 10.**( (AL'JCIO(YTT)-I . )/3.419) 

GO TO C7 

e6 F=1C,»*( ((ALOG10(Y7T)-ALOC10( .4))/1.4367)-^»LCG1C(1.5)) 

67 iU'TP = (G* (1 .-LI )/UF)*F»*1 .25 

S1 = .P4-( (A.LOG10(FFTP )-«,301 )/1.332U) 

TW-TF( 1) 

77 TW=Ti;*.0 1 

I'’ (T2.GT.TS) 00 TO 307 

C=(H’tUf)* (TW-TK1 ) )•* .<?9*ri + ;iC»( 1.-M )••. e*) )*(TV-T’'( 1) ) 

i»='U,TA*(TS-TW)*».75 

i.=Ai'S(C-I>) 

I) (F.CT.25.) GO TO 77 
C1=C/(TJ-TF( 1)) 

L1=D/(T.G-TW) 

Ur = (Cl»D1)/( C1*C1 ) 

AP = ( .1»0i')/(U?,*(TS-TF( 1) )) 

A=Ab+A 
aK CONTI*:ilF 

'Ji =CR/((Tf.-TF(1))«A) 

AT=A+A‘!(> 

WPITL (5,1) 

1 E05XAT (• TS’,5X, 'in*, 5X,*FFF0A',SX,'Z*,7X,’UNF*,7>, •;){>•, 

19/ ,*AM;-' ,9/, ’A' • ,9X, 'AT' ) 

WFTTF (5,275) TS ,TF ( 1 ) ,LFFO), , : , UXL ,0E , ANb, » , AT 



-no- 




t « 



tj 




n 10./,; X ,Kio. ; Flo. D section g 

Z i>i. 1 Cs);.i _ f 

r;' j: = riN 

44 T' :'’:=t'fu7+.oi 

;= (T 0 U 7 +TIN )/. . 

oU=4.13-.0‘)7<5»(TAV-??.) 

=. 314+ .00066* (TAV-3?.) 
vt A0*( 1.-t FFOA) 

'■»=?F/(TfiU I-ll.. ) 

c* = (4.*:‘«)/(r3**2*:'i. n«i!59*m36(- . ) 

V2=Ox//24?tjO . 

; += . 0? 3* ( y ■»/ L 3 ) * ( (;+ * D3 /'j w)**.p*(u«/<w)**. 4 
T'i=T3V 

40 "T=T'r+.0C5 

IF (Tl.GI.TFFf )) 0'' TO 307 SECTION H 

IF (TT.LT.T'V) r,( TO 3C7 
FA :r = HFC<S)*.426*(:F(5)-T'T) 
a6 = 1./VF(S)-1./VC(=,) 

;6=i./vF(s) 

r/=lK6)-TT 

hCUM = (1 .♦(4e.31/FF(;(S))*{Ti {‘>)-TT))*72.9» 

1{( (1./VF(S))*( l./Vi (S)-1./V':<'3))«(KF1«*3)*LNTH)/(r'?*UF1* 

, (Ti (*i)-TT)))**.25 
F;=I; -<'*{"’T-TA V)-.iC0.»r'*(TF(5)-TT) 

L3=iao(L2) 

li (L3.GT.20.) V.C TO 40 
Ur=(HW*KCOND)/(Ht.*!ICOND) 

T0'JTr=2. *(TF(5)-{;F/{UC*A2) ))-TIN 
:.1=A! F.(TCUT-10U71 ) 

TC,UT=TCUT+(TOtJTL»Tf UT)/10. 

IV (n.CT..5) GO TC 44 

01 ={ 1.06 792;- 11)* (GV. ••?)•( (C ..*D3)/UV)**-, 25 
pp=( l,.p»Fk/o2.3 )•( 144 ,/778. )*<1./3.413E 06) 

c ov;;,- ALL THFi-firr' rrrTCTPV'cv 

ETrt= (3.4 13?. C6*( 1000.-PP) )/CA!. 

C rJ. SPLIT 

H:20=(CAli*FfF'^T)/{j.4l3c Of.) SECTION I 

P*V'3 = 1000.-rH20 
WiITt (5 ,444) 

444 F(;PMAT CO',* n*,7X, ‘TIN*, 5X, ‘TOUT*, 4X,*FFF*’,fX,*VW,4X, 

1 ‘PriPO* ,5X, TNi.3* ,3.':, *hHV • ) 

WFITF. (5,4b0) PP.TTF ,T0i:T,FTA,VW,pi^20,FV'H3,Pr 

450 F0R3AT ( F7 . 2 , 2X , F6 . 1 , 2X ,Ff . 1 , 2X ,F7. 5 , 3X ,F6 . 1 , 
12X,F6.1,2(,f6.1»2X,I7.1) 

c COST £;o;.rio*i 

CAPT = !i 1 »A.r-F M!3*( P2-P3 ) 

CApl=CAFT*i9*(IF*((1. + lF)** -4)/(( 1 . .*1® ) * * R9- 1 . ) ) 

‘.AVls((1.47491 10*F6*RP*P9)/(HA*F(F*'’4»H7))*(1./P'^-1./ETA)-CAPT 

«t<ITi (5,451) 

451 Fr5''AT ('O*,* TNTFKIST” ,?x, ’HFATFa S/FT?’ ,2X, ’C/KV H20',2X, 
I’S/'a :.H j’ ,2X, ’POlLtF tFF’,2V,’AUX '■'FF’) 

»MT1 (5,45?) ■!K,F1,R?,h3,PU,HA 
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APPENDIX II STEAM CYCLE DATA 

The computer program discussed in Appendix I requires the 
following information about the steam cycle as inputs: 

(1) The steam condensing temperature, T^. 

(2) The steam cycle thermal efficiency at T 5 . 

(3) The portion of the steam side heat transfer (condensing) 
coefficient in the steam condenser/ammonia boiler, 
which is a function of T^ i.e. the fluid properties. 

Carrying out the modification of the 'Bull Run' plant detailed 
in chapter five, the steam cycle thermal efficiency versus the steam 
condensing temperature is calculated. 

The properties involved in the calculation of the steam side 
heat transfer coefficient and their relationship are shown in 
equation (23). The portion of equation (23) which includes only the 
fluid properties is called 'Bq'- 

Table A2-1 shows values of thermal efficiency, nthc^- '^ 0 ' 

S L6dm 

for the different values of steam condensing temperature used in the 
analysis. 



T$(°F) 


'^^'^steam 


6 o(BTU/(hr-ft- 


.(op)-^^)) 


91.7 


.4578 


3438.6 




no 


.4451 


3642.4 




127 


.4336 


3760.3 




131 


.4308 


3791 .7 




147 


.4194 


3917.2 




162.2 


.4084 


4085.5 




182.2 


.3952 


4270.1 




TABLE A2-1 


VARIATION OF STEAM CYCLE 


EFFICIENCY AND 'Bq' WITH 


THE STEAM 



CONDENSING TEMPERATURE 
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